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Spin susceptibility of superfluid 3 He-B film with specular surfaces is calculated. It is shown 
that, when the magnetic field is applied in a direction perpendiculr to the film, the suseptibility 
is significantly enhanced by the contribution from the surface bound states. No such enhance- 
ment is found for the magnetic field parallel to the film. A simplified model with spatially 
constant order parameter is used to elucidate the magnetic properties of the surface bound 
states. The Majorana nature of the zero energy bound state is also mentioned. 
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Existence of low-energy midgap states in the vicinity 
of the surface and/or interface is a universal feature of 
unconventional superconductors and superfluids. 1-4 ' The 
surface bound states govern the transport properties of 
superconductors and superfluids, since the system always 
communicates with its environment through the surface. 
For example, the zero-bias conductance peak (ZBCP) in 
the tunneling spectrum of high T c superconductors 4-6 ' 
has been ascribed to the zero-energy bound states. ZBCP 
has been also observed in other superconductors, such as 
Sr 2 Ru0 4 , 7) K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br, 8 ) UBei 3 , 9) 
and CeCoIn5, 10 ) and was regarded as evidence of their 
■ unconventionality. 

The surface bound states in p-wave pairing superfluid 
3jj e i,2, 11-13) hayg no ^ been observed until recently, be- 
cause of the lack of an appropriate probe for the neutral 
superfluid. It has been recently reported, 14 ~ 18 ) however, 
that the transverse acoustic impedance Z provides useful 
information on the density of states of the surface bound 
states in superfluid 3 He-B. The surface bound states in 
this system have recently attracted much attention as 
\ Majorana fermion surface states which characterize the 
topological symmetry of the bulk BW state. 19 ~ 24 ) Re- 
cently Chung and Zhang 23 ' suggested that, when the 
spin quantization axis is taken parallel to the surface, 
the low energy behavior of the field operators looks like 
that of Majorana Fermions, as a result, the magnetism 
of surface bound states is Ising like polarized only in the 
direction of surface normal. It follows that the suscepti- 
bility will become anisotropic. They proposed to detect 
this anisotropy by measuring the spin relaxation rate of 
an electron that forms a bubble trapped near the 3 He 
liquid free surface. 

In this letter, we discuss the spin susceptibility of su- 
perfluid 3 Hc-B film. The properties of films with width of 
several coherence lengths are expected to be dominated 
by the surface properties. We show that the susceptibility 
indeed shows anisotropy. When the magnetic field is per- 
pendicular to the film, the susceptibility is significantly 
enhanced. The susceptibility of thin films even exceeds 
the Pauli susceptibility. However, the spin susceptibility 
at T — remains finite and does not diverge as Ising 



spin systems. When the magnetic field is parallel to the 
film, no enhancement of the susceptibility is found. 

Let us briefly review the surface bound states. We con- 
sider a superfluid 3 He-B filling z > domain with a 
specularly reflecting plane surface located at z = 0. The 
Hamiltonian of the system in 4-dimensional Nambu rep- 
resentation is 



«"5 



drdr'¥(r)£(r,r')^(r'), 



where 



#(r) = 



w 



(1) 



(2) 



is the Fermion field operator in the Nambu representa- 
tion. The energy matrix £{r,r') in the presence of mag- 
netic field is a 4 x 4 matrix given by 

f£(V)5(r-r') A(r,r') 



£{r,r') = 



\ At(r',r) 



£(V)<S(r 
6(r-r') 



(3) 



with 7 the gyromagnetic ratio of the 3 He atom and <x 
the transpose of Pauli matrix er. The order parameter 
of the BW state with a surface is given by 



A(r,r') = V e lp - ir - r,) A(z,p) 



p 



A{z,p) 



A (z)(-p x + ip y ) 
A l {z)p z 



Ax{z)p z 
A (z)(p x + ip y ] 



(4) 
(5) 



where p is a unit vector along p. We seek surface bound 
states solving the Bogoliubov equation 



J dr'£{r,r')^{r r ) = E^{r) 



(6) 



under the boundary condition \I> = at z = 0. 

We assume here that An and Ai are constant, for sim- 
plicity. Using the quasi-classical approximation we can 
obtain the bound state wave functions (Nambu ampli- 
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tudes). 15 ' 

Let us first consider the case without magnetic field. 
We can find both positive and negative energy bound 
state for each Fermi momentum (K» , k), where K» is the 
component parallel to the surface and k is the perpendic- 
ular component. The bound state energies are given by 
±Ao sin^ with the polar angle of the Fermi momentum 
with respect to the z axis. 12 ) For positive energy eigen 
value E = Aosin^, the eigen function is given by 



*&i 0) = e iK r r u{k,z) ($+ - e**$_) 
and for negative energy eigen value E = — Ao svctO 
^ { -\r) = e lK r r u {k, z) (e"^$+ + $_) , 



(7) 



(8) 



where (j> is the azimuthal angle of the Fermi momentum 
around the z axis. Here $± are the Nambu amplitudes 
given by 






V-v 



i 
1 



(9) 



It is worth noting here that <I>± is the eigen vector of spin 
operator S z in the Nambu representation 



S z $± =±2 $ ±- 



(10) 



S - l - h 

The z dependence of the eigen functions is included in 

u(k, z) = ue~ KZ sinfcz (11) 

with k = Ai/vp and the normalization constant u is 

determined so that are normalized. 

Since all the eigen functions of £(r, r') form a complete 
set, we can expand the Fermion field operator ^: 

*(r) = £ (7«r„*^(r) + 7^,*^ (r)) + ■ ■ ■ , (12) 
K \\ 

where we have omitted the gapped modes. Since all the 
eigen functions are orthogonal to each other, we obtain 



= 



7k 



J dre- lK r r u {k,z) 

x $ + |* > -e" 
J dre lK r r u (k,z) 
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(13) 



(-e** < $+|* > + < $_|* >) , (14) 



where 



< $+|* >=V>t(r) + iip\(r), 



(15) 
(16) 



< >=(-i)tp l (r) + ip\(r). 

We can show that jK n and 7]^ satisfy the Fermion com- 
mutation relation 



{7*|| , 1K[ } = 0, 



Some caution is necessary about the doubly degener- 
ate zero energy states which happen when Kn = 0. In 



this case, the system is essentially the polar state 2 -* and 
the azimuthal angle <p is an irrelevant quantity. We can 
choose any linear combination of as an eigen func- 
tion of the zero energy state. An example is 

* = (70$- + la^+^j V2u(kp, z)+nonzero energy states. 

(18) 

In this case, 70, 7o are still Fermion operators because 
$± are mutually orthogonal. Another choice is 



r+ ($_ + $+) + r_- ($_ - $+) u(k F , z) 



(19) 

with T+ = (70 + 7 u )/v / 2,r_ = *(7o - ll)/V2. The new 
operators T± have a Majorana property 

r 1 f = r+, rl = r_ (20) 

and obey the commutation relation 

{r + ,r+} = {r_,r_} = i, {r+,r_} = o. (21) 

Now let us consider the effect by magnetic field. To 
obtain the low energy spectrum, we consider matrix ele- 
ments of £ of Eq. (3) between the eigen functions given by 
Eqs. (7) and (8). It is quite interesting that only S z has 



a finite matrix element between and ^ 

" ii 



/(") 

K„ 



Other 

spin components S x , S y have no matrix element at all. 
It implies that the surface bound states respond only to 
the magnetic field in the direction of the surface normal. 
This agrees with the recent suggestion by Chung and 
Zhang. 23 ) When the magnetic field is applied in the z- 
direction, the surface bound state wave function is given 
by 



-K II -til 



(22) 



and the energy is obtained from an eigen value equation 



A sin 9 



\ 




E 



to be 



£7 = ±4/(A sin^) 2 + 



(23) 



(24) 



In the If || — ► limit, the wave function (22) is reduced to 
$_ for the positive energy and to $ + for the nega- 
tive energy Let us consider the spin susceptibility 
of the ground state. Contribution from the occupied neg- 
ative energy bound states to the ground state energy is 
given by 



Kn 



"(A (1 ,in^-> + ( ^ 



(25) 
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2 V2A 



(26) 



where the factor 1/2 in Eq. (25) comes from the prefactor 
in Eq. (1). We obtain the susceptibility 



Xz 



d 2 ~f 2 k 2 F 
g^Eo(H) 



(27) 



16ttAo 

which is as large as the normal state susceptibility Xn 
multiplied by the width l/n — vp/Ai of the bound 
states. The susceptibility of the surface BW state at 
T = is large but finite, while in the polar state the 
susceptibility will diverge because A = 0. The mag- 
netism of the polar state is just Ising spin like because 

the bound state energy under magnetic field splits into 
+ 1 H 

- 1 - 2 ' 




Fig. 1. Self consistent order parameter at T = 0.2T C of 3 He-B 
film with width L = 5, 6, 10, 15, 20§. The coherence length £ is 
defined by § = vp/irT c , which is equal to 0.155/^m at SVP. 



Now we turn to the susceptibility of films. To calculate 
the susceptibility, we use the quasi-classical Green func- 
tion method developed for dealing with the boundary 
problems. 12, 25 ~ 27 ) Within the quasi-classical Green func- 
tion theory, however, it is not straightforward to calculate 
the linear response in finite systems because the quasi- 
classical Eilenberger equation needs an initial condition 
at some point. We instead start from the linear response 
formula for the Gor'kov Green function and rewrite the 
result using the quasi-classical Green functions and the 
evolution operator. 26,27 ) When we apply the formulation 
to the case with the constant order parameter An, Ai , we 
obtain the susceptibility ( per unit volume ) Xzz (0) at the 
surface 

r ir/2 

(28) 



Xzz(0) - Xn 



XN 



d6sm9n;T V, 

n>0 



A 2 cos 2 e 



K + A 2 sin 2 6)^1 + A 2 si 



sin 2 9 + A 2 cos 2 I 



(29) 



where xn is the normal state Pauli susceptibility and w n 
is the Matsubara frequency. 

In the presence of surfaces, the order parameter is 



modified cdas can be seen in Fig. 1. 12 > 13 ) Near the spec- 
ular surface, the perpendicular component A\(z) is sup- 
pressed to zero, while the parallel component An(z) is 
somewhat enhanced to compensate the pairing energy. 
In this report, we show the results of numerical calcula- 
tions using the self-consistent order parameter of Fig. 1. 
The details of the calculation shall be reported elsewhere. 




Fig. 2. Local distribution of \ zz in 3 He-B film at a temperature 
T = 0.2T C . \zz is scaled by the normal state value xn ■ Fermi 
liquid correction by is taken into account. Film widths are 
the same as in Fig. 1. 
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Fig. 3. Local distribution of in ^He-B film at a temperature 
T = 0.2T C . Xxx is scaled by the normal state value xjv- Fermi 
liquid correction by Ffi is taken into account. Film widths are 
the same as in Fig. 1. Note that the vertical scale is different 
from Fig. 2 



In Fig. 2, we show the local distribution of Xzz( z ) a t 
a low temperature T = 0.2T C for several width films. 
The Fermi liquid correction by Fq is taken into account. 
The vertical axis is the local susceptibility normalized 
by the normal state susceptibility. The enhancement of 
the susceptibility can be clearly seen at the end surfaces. 
The bottom value is almost equal to the bulk B-phase 
susceptibility. In sufficiently thin films, the overall en- 
hancement is found rather than the surface enhancement. 
This is because the bound state wave functions extend 
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over the entire width of the film. In constract to Xzz, the 
susceptibility x.xx for the magnetic field parallel to the 
surface is not enhanced at the surfaces as can be seen 
in Fig. 3. The surface value of x.xx is even smaller than 
the bulk susceptibility and is nearly equal to that of the 
planar state with the Fermi liquid correction. These re- 
sults clearly demonstrate that the surface bound states 
respond to the magnetic field only in the direction of the 
surface normal. 

Finally we show in Fig. 4 the susceptibility Xzz and 
Xxx averaged over the film width. We can find that \zz 
even exceeds the normal state Pauli susceptibility for suf- 
ficiently thin films. In thicker films, Xzz is still larger than 
the B-phase bulk value, while Xxx remains smaller. The 
anisotropy is large enough to be observed. 




Fig. 4. Film width dependence of \ zz and Xxx averaged over the 
film width at T = 0.2T C . 



In conclusion, we have shown that the spin suscepti- 
bility of superfluid 3 He-B film has a strong anisotropy 
caused by the surface bound states. The anisotropy is 
sufficiently large to observe experimentally. We have con- 
sidered, however, films with specular end surfaces. In 
actual films on the substrate, the surface scattering by 
the substrate will be diffusive. One of the methods to 
avoid the diffuse scattering is to coat the substrate by 
He layer. In fact, recent experiments of the transverse 
acoustic impedance 17 ' 18 ) showed that the specularity of 
the surface is considerably enhanced by the coating. On 
the other hand, the susceptibility of the film with diffu- 
sive surfaces is itself of interest. The density of state at 



zero energy is known 11-13, 24 ) to be increased by the diffu- 
sive scattering, which might lead to further enhancement 
of the susceptibility. 
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